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ABSTRACT 

We report the resuhs of Keck i'-band non-redundant aperture masking of HR 8799, a system with 
four confirmed planetary mass companions at projected orbital separations of 14 to 68 AU. We use 
these observations to place constraints on the presence of planets and brown dwarfs at projected 
orbital separations inside of 10 AU — separations out of reach to more conventional direct imaging 
methods. No companions were detected at better than 99% confidence between orbital separations 
of 0.8 to 10 AU. Assuming an age of 30 Myr and adopting the Baraffe models, we place upper limits 
to planetary mass companions of 80, 60, and 11 Mjup at projected orbital separations of 0.8, 1, and 
3-10 AU respectively. Our constraints on massive companions to HR 8799 will help clarify ongoing 
studies of the orbital stability of this multi-planet system, and may illuminate future work dedicated 
to understanding the dust-free hole interior to ~6 AU. 

Subject headings: instrumentation: adaptive optics — instrumentation: interferometers — planets and 
satellites: detection — stars: individual (HR 8799, HD 217165) — techniques: high 
angular resolution 



1. INTRODUCTION 

HR 8799, a 20-150 Myr (|Mo6r et al.l 120061) A5V star 
hosting several planets, presents a challenge for for- 
mation models of massive exoplanets. Until recently, 
this system was known to have three ~5-7 Mj^p ex- 
oplanets (assuming an age of 30 Myr), clearly identi- 
fied with ada ptive optics ( "AO" ) at separations of 24- 
68 AU ((Marois et al.ll2008D . The recent id entification of 
a fou rth ~7 Mjup companion at 14 AU dM arois et aLl 
120101) reinforces that this system contains an intrigu- 
ing architecture. Maintaining this system in a dy- 
namically stable, non-resonant configuration for the 
10^-10* year age of the system has proved to valid 
only for models poss essing a restricted range of physi- 
cal parameters, (e.g. iGozdziewski fc Migaszewskil 120091 : 
iFabrvckv fc Murrav-Clavll2010D . Furthermore, forming 
massive planets such as these at such wide projected 
orbital separations is problemat ic for the core-accretion 
model of planet formation (e.g. IDodson-Robinson et al.l 
[20091) . 

The recent direct imaging observations dedicated to 
HR 8799 are most sensitive to planetary mass compan- 
ions at projected orbital separations >10 AU. However, 
obtaining a census of any planets within 10 AU is essen- 
tial for a full und erstanding of the origin and stability 
of the system (e.g. lReidemeister et al.|[2009l) . Several au- 
thors have addressed the possibility that the placement of 
wide-separation exoplanets is caused by outward scatter- 
ing f rom inner m assive perturbers or a planetesimal disk 
(e.g. lVeras fcA rmitagc 2004; Vcras ct al. 2009). The ad- 
vantage of such models rests on the formation of the bod- 
ies just beyond the location of the ice line at a few Myr 
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where the core accretion timescales are much shorter. In 
a similar vein, several works have analyzed the degree to 
which the placement of wide separation planets such as 
those in HR 8799 are caused by mutual outward sca ttcr- 
ing from t he existing planets (e.g.lChatteriee et al.l2008.,) . 
Similarlv. iScharf fc Menoul ()2009() justify the existence of 
long period exoplanets as a natural by product of dynam- 
ical relaxation through planet-planet scattering. 

Further, several multi-planet exoplanetary systems ap- 
pear to be "dynamically p acked" (Barnes fc Ravmondl 
[20041 : [Ravmond et al.|[2009D : i.e., with the placement of 
their mutual orbits lying close to a boundary for dynam- 
ical stability. Hence, in systems with multiple compan- 
ions, it is imperative to assess the full orbital distribution 
of planetary mass companions and test such hypotheses. 

The planets orbiting HR 8799 discovered to date have 
been imaged using ang ular differential ima ging ("ADI" — 
Lafreniere et al.l [20071: [Marois et aI|[2006tlLeconte et al.l 



20101) . To obtain high contrast images of this system at 



even smaller separations, we have used non-redundant 
aperture masking, hereafter "NRM", using the NIRC2 
infrared camera at the W. M. Keck Observatory. NRM 
is the only technique that can produce significant con- 
trast levels at such small inner working angles (< 300 — 
500 mas) among existing 10m class telescopes. We are 
able to place upper limits on additional massive com- 
panions to the HR 8799 system that would otherwise 
have been veiled by the telescope diffraction or the as- 
sociated uncorrecte d quasi-static speckle noise (see e.g. 
[Hinklev et aI1[200g)^ using more conventional dire ct imag- 
ing techniques (|Oppenheime~ fc HinklevI l2009l) . Our 
technique, which we describe in greater detail in Sec- 
tions|2]and|31 bolsters ongoing ground-base d NRM efforts 
fireland fc Kraus .2008: .Kraus et al,, ,20081: [Bernat et all 
[2010D as well as serving as a precursor for the NRM 
efforts to be used with the James Webb Space Tele- 
scope, hereafter "J WST" ([Sivaramakrishnan et al.[|2"Q10l : 
[Dovon et ai]|2010l) . 
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2. METHODS & OBSERVATIONS 
2.1. Non Redundant Aperture Masking 

NRM (jBaldwin et all [l986l: iReadhead efall [l988[) 

achieves contrast ratios of '^10^ - 10'^ at very small in- 
ner working angles, usually within ^^X/D — A\/D 
20-300 mas for Keck L'-band im aging) . Applicatio ns of 
this technique fe.g-.lriltliill et al.ll200aiLlovd et al.ll2006l: 



A summary of the observations is presented in Table [T] 



Monnier et alll2007t IWoodruff et al.ll2008HBernat et alj 
20101) use AO along with an opaque mask containing sev- 



eral holes, constructed such that the baseline between 
any two holes samples a unique spatial frequency in the 
pupil plane. A key feature of this technique is its ability 
to sample the " closure phase" quanti t y between any thre e 
baselines (e.g. [Baldwin et all Il986l: iHaniff et al.1il987D . 
While the effective transmission of the telescope is re- 
duced by ^90-95% due to the mask, the effectiveness of 
NRM is limited primarily by the ability to calibrate the 
point spread function. Hence the importance of measur- 
ing closure phases using carefully chosen calibrator stars 
exceeds the need for gathering additional flux. Assum- 
ing the aperture holes are small relative to the charac- 
teristic size of any atmospheric turbulence, and in the 
regime where the AO system is providing overall phase 
stability, the closure phase quantity is largely unaffected 
by atmospheric turbulence. Under good conditions, this 
technique can achieve L'-band contrasts of 10^-10^. This 
technique is particularly well suited for stu dies of stellar 
multiplicit y and the bro wn dwarf desert (e.g. lKraus et all 
[2008: Ireland fc KrausI 12008) . or detecting giant planets 
orbiting young stars (< 50-100 Myr). 

2.2. Observations 

We imaged the HR 8799 system using NRM operating 
at L'-band (3.76 ^m) on the nights of UT 2009 Aug 5 
and 6 using AO and the NIRC2 infrared camera on the 
Keck 2 telescope at Mauna Kea. While no seeing mea- 
surements were made, non-NRM observations obtained 
each night verified that the AO system was able to return 
apparently diffraction-limited images at K-hanA. An ob- 
serving sequence consisted of observing a star in both the 
upper left and lower right quadrants of the NIRC2 detec- 
tor with the Keck nine hole aperture mask (see Figure 
1, left panel). At each detector position, we obtained 15 
images with 20 s exposures. Nine and ten such 30-image 
sequences of HR 8799 were obtained on the first and sec- 
ond nights, respectively. Observations of a calibrator star 
were obtained using a similar 30-image sequence in be- 
tween each of the HR 8799 target sequences, with the cal- 
ibration star interferogram placed in the identical man- 
ner on the array. A different calibrator star was observed 
in between each HR 8799 sequence. The HR 8799 and 
calibrator star observations, plus observing overheads, 
required ~4 hours on each night. We selected these cal- 
ibrator stars from the compilatio n of stars with stabl e 
radial velocity measurements by INidever et al.l (j2002D . 
who found them to have RMS variations of <0.1 km s~^ 
over an interval of ^10 years; this criterion was meant to 
select against binary systems that would not be suitable 
for calibration. All calibrators were selected to be near 
HR 8799 on the sky {p < 15°) and to have similar optical 
V magnitudes (to obtain similar AO correction) and sim- 
ilar or brighter near infrared K magnitudes (to match or 
exceed the signal-to-noise in the HR 8799 masking data). 




-2-10 1 2 -2-10 1 2 
arcsec Spatial Frequency (10^ rad'^) 

Fig. 1. — Left: A Keck L'-band interferogram of the star HR 
8799 produced with the NIRC2 nine hole non redundant aperture 
mask. Right: The power spectrum corresponding to the image on 
the left. 



3. ANALYSIS 

To use the closure phase quantity to search for plan- 
et ax^;_jnBSS_compajQions, w e follow the analy sis outlined 
in IKraus et all pOOl and [Ireland fc Krausi (pOOl . We 
briefly summarize the procedure here. The data are 
initially flatfielded, sky subtracted, aligned, and cor- 
rected for cosmic rays. The bispectrum, the complex 
triple product of visibilities defined by the three baselines 
formed from any three subapertures, is then calculated. 
The phase of this c omplex quan t ity is the closure phase. 

As discussed in IKraus et al.l (|2008| ). the calibration 
procedure centers around calculating the closure phase 
for calibrator stars. The calibrated object closure-phase 
is found by subtracting a weighted average of the closure- 
phase for the calibrator stars. For the analysis in this pa- 
per, which is motivated by the search for point sources, 
the squared visibilities were not used as they were nois- 
ier than the closure-phases. Eight separate calibrators 
were used, as detailed in Table [TJ Each calibrator was 
calibrated against all other calibrators to search for re- 
solved calibrators. HD 217165 was found to be a binary 
star of separation 125 mas and 2.6 magnitudes contrast 
during this process, and was removed as a calibrator. 
Each target data set was then calibrated with each cal- 
ibrator observed on the same night. The Root Mean 
Square (RMS) calibrated closure-phase was found for 
each of these target-calibrator pairs, and all calibrations 
that resulted in an RMS closure phase more than 1.5 
times the minimum for each target over all calibrators 
were assigned a weight of zero. In practice, this meant 
that each target data set was calibrated by an average 
of the two or three calibrator data sets obtained closest 
in time. The causes of the deviations from zero closure- 
phases (of <2 degrees RMS) in the calibrator stars and 
the time-dependence of calibration fidelity are not fully 
understood, but we have evidence that dispersion in the 

^ Angular separation between the calibrator stars and the HR 
8799 target, measured in degrees. 

^ Total exposure time: number of exposures multiplied by the 
exposure time for each image for UT 2009 Aug 5. 

3 Total exposure time for UT 2009 Aug 6. 
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TABLE 1 

Observing log for. HR 8799 and Calibrators. 



Object K Separation (deg) ^ RA (J2000) Dec. (J2000) texp (s) ^ texp (s) ^ Type Comment 



HR 8799 


5.24 




23 07 28.7 


+21 08 03.3 


270x20 


300x20 


target 


HD 217813 


5.15 


1.12 


23 03 05.0 


+20 55 06.9 


30x20 


60x20 


calibrator 


HD 216625 


5.73 


3.56 


22 54 07.4 


+19 53 31.4 


30x20 


30x20 


calibrator 


HD 219172 


6.09 


5.98 


23 13 48.2 


+15 22 03.5 


30x20 


30x20 


calibrator 


HD 218133 


5.68 


6.70 


23 05 31.2 


+14 27 06.0 


30x20 


30x20 


calibrator 


HD 217165 


6.19 


11.53 


22 58 29.9 


+09 49 31.9 


60x20 


30x20 


calibrator binary (not 


HD 222033 


5.79 


12.08 


23 37 06.6 


+30 40 40.9 


30x20 


60x20 


calibrator 


HD 221830 


5.30 


12.11 


23 35 28.9 


+31 01 01.8 




30x20 


calibrator 


HD 210460 


4.47 


14.37 


22 10 19.0 


+19 36 58.8 


30x20 


30x20 


calibrator 



L'-band due to water-vapor is one cause. The lack of 
perfect cahbration is stih the dominant source of closure- 
phase noise in this analysis, however the final contrast 
limits are within 1 magnitude of being photon-limited. 

4. RESULTS 

In this section, we search for companions in the the HR 
8799 dataset by identifying any non-zero closure phases, 
to be elaborated below. No companions were detected 
interior to 14 AU, the location of HR 8799e. We thus 
use these data to place limits on companion masses with 
orbital radius. 

4.1. Detection Limits 

To find the contrast limits of our data, we first added 
a calibration error in quadrature to uncertainties esti- 
mated directly from the scatter in the data so that the 
chi-squared for a null model with no companion was 
equal to the number of measured closure-phases. We 
then scaled the chi-squared value by the ratio of the to- 
tal number of linearly independent closure phases to the 
total number of closure phases (28/84) for the purpose 
of standard least squares fitting, which assumes inde- 
pendent measured quantities. We have found this tech- 
nique to be nearly equivalent but much simpler than 
the detai led estimation of covariance matrices as de- 
tailed in iKraus et all ()2008[ ). We then simulated 10'* 
data sets with closure-phase standard deviations that 
matched these modified uncertainties. For each simu- 
lation, we found the best contrast at each value of sep- 
aration (searching over position angle) and assigned the 
99% confidence limit to the contrast where 99% of the 
simulations had no best fit companion above this limit. 
When fitting to the actual data using a grid search, we 
also found no tentative solutions above this limit, there- 
fore we are reporting a null result. 

Figure [5] shows our L'-band detection limits achieved 
with the two nights of Keck NRM data. In Figure [H 
we show the brightnesses and positions for the four 
planetary rn ass companions discussed in iMarois et al.l 
(pool mmh . The HR 8799 "d" and "e" companions 
are within the angular range probed by our NRM obser- 
vations. Both components are approximately AL'~9.3 
mag fainter than the star, while our observations achieve 
AL'~7.78 and AL'~7.74 at these orbital separations and 
are unable to detect these companions. However, we 
are able to probe much closer in to the host star than 
these separations. We achieve our best contrast of 7.99 
mags at the L'-band diffraction limit of the Keck tele- 
scopes, corresponding to a 3 AU projected orbital sep- 
aration for HR 8799. This orbital separation is compa- 



rable to the 2-4 AU ice lin e boundary for a mid-A star 
([Kennedy &: Kenyonll2008[ ) in the 3-5 Myr timeframe for 
formation of gaseous giant planets. 

To convert the contrast limits to constraints o n com- 
panion masses, we used the models of Baraffe et al.l 
( 2003 ) along with the Hipparcos distance of 39.4 pc 
(|van Leeuwen|[2007l ) ■ These models were chosen to main- 
tain consistency with IMarois et al.l (|2008l . I2010[ ) as well 
as to provide ease of comparison with other studies. We 
assume an age of 30 Myr, which is consiste nt with the 
20-150 Myr age assigned to this system by I Moor et all 
(|2006D as weU as th e two age regimes (3 01 ^j] and GOtlg^ 
Myr) discussed in IMarois et all (|2010( ). (iMoor et al.l 
[2001 assigned an age of 20-150 Myr to the HR 8799 
system. In addition, there is evidence that the sys- 
tem may be a member of the Colu mba moving group , 
and hence has an age of 30 Myr ([Torres et al.l l2008l 
R. Doyon, B. Zuck erman, private communication 2010). 
iMova et all ([2010D . however, suggest the system may be 
as old as 1 Gyr. For consistency withTMarois et al. (200l, 
I2010( ). we adopt an age of 30 Myr. 

Table[5]summarizes the sensitivity limits to companion 
masses for an assumed age of 30 Myr. At a separation 
of 3 AU, we place an upper limit of 11 Mjup to any 
companion, with a comparable sensitivity outwards to 
24 AU. Interior to 3 AU, we can place upper limits to 
additional companions of 13 Mjup and 60 Mjup at 2 AU 
and 1 AU, respectively. Our measurements are also able 
to probe to sub-AU levels in the system placing an upper 
limit mass of --80 Mjup at 0.8 AU. For an age of 100 Myr, 
the masses of the currently known exoplanets would be 
estimated at 9-14 Mjup, and our sensitivity at 1, 2, and 
3 AU shift to 80, 30, and 22 Mjup, respectively. 

5. CONSTRAINTS ON INNER MATERIAL: NO MASSIVE 
PERTURBERS 

As discussed in the Introduction, planetary scattering 
is one possible model to explain the presence of four mas- 
sive planets at large orbital radii in the HR 8799 system. 
An object interior to, and more massive than, the four 
5-7 Mjup planets would likely be nec essary to scatter 
the o bjects to their current placement ([Chatteriee et alj 
|2008( ). Our observations place the most stringent con- 
straints to date on any massive perturber interior to the 
orbit of HR 8799 "e." The mass of any perturber must 
be less than ~11 Mjup between 3 and 10 AU. It is pos- 
sible, of course, that a perturber more massive than our 
limit was in fact present at 3 AU, but has since moved 
significantly inward of 3 AU or even was scattered into 
the star. 

The dynamical stability of this four planet sys- 
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Fig. 2. — Detection limits for the two combined Keck nights of non-redundant aperture masking data on the HR 8799 system. The 
green line shows the minimum brightness required for a 99% confidence detection in our HR 8799 data. The dotted orange lines show 
selected corresponding mass limits as described in BarafTc et al. (2003) assuming an age of 30 Myr. Also shown is the theoretical diffraction 
limit at L'-band on a 10m telescope (vertical dotted line), as well as the L'-band brightnes ses of the four known companions to HR 8799 
IjMarois et al.|[2003 . 1201CI '1 . The blue regions represent the debris disk structures defined by ISu et al] 1)20091 ') . 



tern has been discussed a t length i n other works 
(iGozdziewski fc Migaszewskil [2009: Moro-Marti'n et alJ 
120101: IFabrvckv fc Murrav-Clavl I20m) . Stabihtv argu- 
ments almost certainly require that the planets are 
both in a resonant configuration and all have masses 
<10Mjup. Indeed, iMarois et all ()2010[ ) point out that 
only a small fraction (7 out of 10^) of trial dynamical 
simulations of the system remain stable over '^160 Myr. 
In addition, the stud y of th e HR 8799 debris disk car- 
ried out by ISu et alJ (j2009() . indicates the presence of 



an inner warm (~150 K) debris belt (see Figure [2]). The 
outer boundary of this disk is most l i kely sc ulpted by the 
"e" component at 14 AU. ISu et al.l (|2009l ) suggest that 
few, if any, dust grains exist interior to 6 AU, and pos- 
tulate that additional planets may be responsible for the 
clearing. Our measurements show no evidence for com- 
panions more massive than '^12 Mjup at 3-6 AU, but 
they do not rule out the possibility that an object less 
massive than ~12 Mjup was responsible for this inner 
clearing. Indeed, even relatively low mass objects (0.1-3 
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TABLE 2 

NRM Detection Limits and Mass Upper Limits to Additional 
HR 8799 COMPANIONS 



Sep. (AU) 


Ang. Sep. (mas) 


AL' (mag) 


Mass Limit (Mjup) ^ 


0.8 


20 


5.17 


80 


1 


25 


5.63 


60 


2 


50 


7.48 


13 


3 


75 


7.99 


11 


3 - 10 


75 - 750 


7.99 - 7.63 


11 - 12.5 



In the more distant future, NRM wih play a 
Tominent role for the planet-finding efforts of JWST 
Beichman et al.l I2010D , significantly improving on the 
sensitivity presented in our study. The Tunable Filter 
Instrument camera which is integrated with the Fine 



M.Tup) can shepherd debr is disk belts (e.g. iQuillen et al.l 
[200llChiang et al.l 120091) . 

5.1. Future Work 

Achieving very small inner working angles (10-250 
mas) corresponding to projected orbital separations of 
a few AU will be in crucial for the future ground and 
space-based characterization of the circumstellar regions 
around young stars. Reaching appreciable contrast at 
such small inner working angles will allow significant 
sensitivity to objects orbiting much closer to the host 
star. Perhaps more importantly, the small inner working 
angles will allow observers to probe solar system scales 
of more distant targets, thereby increasing the num- 
ber of available target stars, includin g those youiig star s 
in the nearest star forming regions iKraus et al.l 120081 ). 
Though traditional high c ontrast techniques such as dual 
imaging polarimetry (e.g. iQppenheimer et al.ll2008D can 
achi eve significant contra sts of ^^8-10 mag within 500 mas 
(e.g. lHinklev et al.|[2009l ). the success of these techniques 
largely depends on a high polarization signal exhibited 
by the planets. However, NRM will have superior sen- 
sitivity at small inner working angles (20-300 mas) and 
will play a major role for ground based high contrast 
imaging of exoplancts and imaging clo se binaries (e.g. 
ILlovd et al.][2006 : Hinklev et al.il2C)llaD . Recent resuhs 
suggest ~7.5 magnitudes of K-h&nd contrast at Keck has 
been achieved for K = 8 star in three hours of observing 
(A. Kraus, private communication, 2010). 

The science return from NRM will be enhanced when 
the imaging science camera has multi-wavelength capa- 
bilities such as an integral field spectrograph (IFS). This 
has already been achieved with the Pro ject 1640 IFS 
and co ronagraph at Palomar Observatory (jHinklev et al.l 
l2011bl Zimmerman et al. 2011, in prep), and will be 
integrated with t he Gemini Planet Ima ger co ronagraph 
(jMacintosh et al.l 120081 : ISoummer et al.ll2009| ). Further- 
more, aperture masking will continue to be an efficient 
tool for characterizing stellar multiplicity. As several 
authors have suggested, more massive stellar systems 
may be more efficient at producing faint stellar compan - 
ions ( !Dodson-Robinson et al.]|200a iKratter et al.ll2010D . 
Aperture masking will play a crucial role to supplement 
ongo i ng surveys of A star mult iplicity (e.g. lHinklev et all 
[2010HZmI^erman et al1[20Tol ). 



Mass limits calculated using the IBaraffe et all 1200^ ') models 
assuming an age of 30 Myr. 

Guidance System on J\V ST fe.g. lSivaramakrishnan et all 
l2010t iDovon et al.ll2010[ ) will incorporate non-redundant 
aperture masking capabilities. Recent simulations indi- 
cate a limiting delta magnitude of ^9.5 for integration 
times of a few hundred seconds in the JWST filters cor- 
responding to the L' and Af-bands. These limits will 
provide sensitivity to objects of a few Jupiter masses for 
syste ms at ages of 10-100 Myr ([Sivaramakrishnan et al.l 
l2OT0h . 

6. CONCLUSIONS 

The results presented in this work are the first con- 
straints on the masses of additional companions to the 
HR 8799 system at projected orbital separations inside 
of 10 AU. We can rule out the presence of a very mas- 
sive perturber that has ejected the four imaged planets 
out beyond their presumed formation region near the ice 
line. Further, for an assumed age of 30 Myr, we place 
upper limits to planetary mass companions of 80, 60, and 
11 Mjup at projected orbital separations of 0.8, 1, and 
3 AU, respectively. These findings will be essential for 
future works to begin clarifying the complicated dynam- 
ical interplay of the four known objects in the system, 
but also furthering our understanding the complicated 
interaction between these companions and the debris disk 
belts in the system. These results may also illuminate fu- 
ture work dedicated to understanding the dust-free hole 
interior to ~6 AU. 
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